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O ' Abstract 

In this paper we take the first step towards a complete next-to-next-to-leading logarithmic 

3 ! 

(NNLL) calculation of the inclusive decay rate for B — > X s j. We consider the virtual correc- 
tions of order a^nf to the matrix elements of the operators 0\, O2 and 0$ and evaluate the real 

<N 

J> ' and virtual contributions to O7. These corrections are expected to be numerically important. We 

t-H ! 

observe a strong cancelation between the contributions from the current-current operators and O7 

O , 

and obtain, after applying naive non-abelianization, a reduction of the branching ratio of 3.9% (for 



/i = 3.0 GeV) and an increase of 3.4% (for /i = 9.6 GeV). 
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I. INTRODUCTION 



Currently, measurements of the inclusive branching ratio BR(B — > X s ^) are provided by 
CLEO [1] (Cornell), by the B factory Belle [2] (KEK), by ALEPH [3] (CERN) and by the 
preliminary BABAR [4, 5] (SLAC) results, leading to a world average of [6] 

BR(5 -> X s7 ) cxp = (3.34 ± 0.38) x 10" 4 . (1) 

This experimental average is in good agreement with the theoretical prediction based on 
the Standard Model (SM) including next-to-leading logarithmic (NLL) QCD corrections 
supplemented by certain classes of leading order electroweak terms [7-10]. For a recent 
status report on inclusive rare B decays and a complete list of references on NLL calculations 
of BR(B — > X s 7) the reader is referred to [11]. In earlier analyses [8, 12-15], the ratio 
m c /rrib, which enters the calculation of the decay width T(B — > X s j) for the first time 
at the NLL level, was tacitly interpreted to be the ratio of the pole quark masses. Using 
m c /m b = 0.29 ±0.02, one obtains BR(B -> X S 7) SM = (3.35 ±0.30) x 10~ 4 , where the errors 
due to the uncertainties in the various input parameters and the estimated uncertainties due 
to the left-over renormalization scale dependence were added in quadrature. More recently, 
Gambino and Misiak [16] pointed out that the branching ratio rises to BR(B — > A s 7)sm — 
(3.73±0.30)xl0~ 4 [16] (see also [17]), if one interprets m c /m b tobem c (/i)/m 6 = 0.22 ±0.04, 
where m c (/x) is the charm quark mass in the MS-scheme, evaluated at a scale \i in the range 
m c < (j, < m b , and m b is the bottom quark IS* mass. 

Despite the current theoretical dispersion on the branching ratio, the agreement between 
the present experimental results and the SM is quite impressive and this has been used to 
derive model independent bounds on the Wilson coefficients C 7 {m w ) and C 8 (mw) (see, for 
example, Ref. [18]). 

Formally, the approximately 11% discrepancy in the branching ratio, stemming from the 
two different schemes for m c /m b , is a NNLL effect. As the measurements of BR(B — > X s ^) 
will become much more precise in the near future, it will become mandatory to systematically 
extend the theoretical predictions to NNLL precision, in order to fully exploit this process 
in the search for new physics. 

To illustrate the complexity of such a calculation, we briefly explain the theoretical frame- 
work. Usually, one works in the effective field theory formalism of the SM, where the W 
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boson and heavier degrees of freedom are integrated out. This results in an effective Hamil- 
tonian in which operators up to dimension six are retained. Adopting the operator definition 
of [12], the relevant Hamiltonian to describe the processes b — > 57, b — > sg and b — > s^g 
reads 

AC 8 

#eff = ^A^C^O^), (2) 

^ 2 i=i 

where Gf is the Fermi coupling constant, X t = V ts *V t b (with V^- being elements of the 
Cabibbo-Kobayashi-Maskawa matrix) and Ci(fi) are the Wilson coefficient functions evalu- 
ated at the scale \i. For practical reasons it is more convenient to use instead of the original 
functions Ci(fx) certain linear combinations, the so-called "effective Wilson coefficients" 
Cf(fx) introduced in [12, 19]: 

Cf(jj) = GO*), (i = l,...,6), 

6 

i=i 

6 

Cf(v) = C 8 (//) + J>C^) , (3) 

i=l 

where yi and Zi are defined in such a way that the leading order matrix elements (s7|Oj|6) 
and (sg\Oi\b) (i — 1, . . . , 6) are absorbed in the leading order terms of C 7 cff (/i) and C 8 efr (yu). 
The explicit values of { Vi } and { Zi }, y = (0, 0, -±, -|, -f , -f ), z = (0, 0, 1, -±, 20, -f ) 
were obtained in Ref. [12] in the MS scheme using fully anticommuting 75 which is also 
adopted in the present paper. 

The operators relevant for our calculation read 

01 = (s Ll ,T a c L ) (c Ll »T a b L ) , 

02 = (sl1^c l ) (c L Yb L ) , 

4 = (5 l7 „T°&l) E(^ T ^) > 
<? 

0s = ii^ ( 5 ^ T ° 6 *) G ^ • ( 4 ) 

Here e = v^ro^ and g s = \[Akol s denote the electromagnetic and strong coupling con- 
stants, respectively. Furthermore, and are the corresponding field strength tensors 
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and L = (1 — 7s)/2 and i? = (l + 7 5 )/2 stand for left- and right-handed projection operators. 
The factor m fe (/i) in the definition of 7 and 8 denotes the bottom mass in the MS scheme. 

For a complete NNLL calculation in this framework, the evaluation of three parts is 
necessary: (1) the computation of the matching coefficients to order a 2 which requires 
a three-loop calculation; (2) the evaluation of the anomalous dimension matrix to order 

where four-loop diagrams are involved; and (3) the calculation of the order a 2 QCD 
corrections to the matrix elements (s^\Oi(ii)\b) (// is of order m h ) which, depending on the 
operator, is either a two- or three-loop calculation. 

The relatively large dependence of the NLL prediction for BR(_B — > X s 7) S m on the scheme 
for m c /mb illustrates that NNLL effects, in particular those related to step (3), can be rather 
large. 

At this point we should stress, that the issue related to the definition of m c /m h serves us 
as a motivation to initiate a NNLL calculation for BR(B — > X^). In the present paper we 
are working out a class of NNLL corrections (to be specified below) to step (3), which is not 
related to the m c /mb issue. However, in many other processes it is known that the kind of 
terms considered in this paper are the source of very important higher order corrections. 

In this paper we consider those corrections of order a 2 to the matrix elements for B — > X s 7 
associated with the operators 0±, 2 , 7 and 0$ which involve a closed fermion loop. It 
is needless to say, that at the same time also the matching coefficients and the anomalous 
dimension matrix should be improved accordingly. Motivated by the fact that the NLL 
corrections to the matrix elements were numerically more important than the improvements 
in the Wilson coefficients, we assume for the time being that this could also be the case 
at the NNLL level. Therefore, we only concentrate on NNLL corrections to the matrix 
elements. In principle also the contributions from the operators 0{ {% — 3, . . . , 6) should be 
considered. However, as the corresponding Wilson coefficients are small, we neglect these 
contributions. Furthermore, we also neglect the NNLL bremsstrahlung corrections to the 
interferences (d,Oi), (O u 2 ), (0 2 ,0 2 ), (O u 7 ), (Oi,0 8 ), (0 2 ,0 7 ), (0 2 ,0 8 ), (0 7 ,0 8 ) and 
(Os,Og), since these terms are infrared finite for vanishing gluon energy and numerically 
relatively small at the NLL level [20]. 

The fermionic corrections we are interested in are essentially generated by inserting a 
one-loop fermion bubble into the gluon propagator of the lower order Feynman diagrams. 
For the numerical evaluation we will assume that rif — 5 massless fermions are present in 
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the fermion loop. 

Once the corrections of 0(a 2 s nf) are available, it is suggestive to use the hypothesis of 
naive non-abelianization (NNA) [21] in order to estimate the complete corrections of order 
a 2 . This is based on the observation that the lowest coefficient of the QCD (3 function, 
fa = 11 — 2nf/3, is quite large and thus it is expected that the replacement of rif by 
— 3/9o/2 may lead to a good approximation of the full order a 2 corrections. There are many 
physical observables, where NNA provides an excellent approximation to the full two-loop 
corrections [22, 23] like the inclusive cross section e + e _ — > hadrons, the hadronic r decay 
or the two-loop relation between the MS and pole quark mass. In particular, we want to 
mention the semileptonic decay T(b — > clvi) where the deviation of the ce 2 /3o terms from the 
complete a 2 result [24] is less than 20%. We also note that the O(a 2 j3o) corrections to the 
photon energy spectrum in B — > X s ^ (away from the endpoint) were calculated in Ref. [25]. 

Our presentation is organized as follows: in Section II we discuss the virtual corrections 
associated with 12 and compute in Section III both the real and virtual corrections to 
Of. The virtual corrections to 8 are considered in Section IV. In Section V we combine 
our findings with the existing NLL results and perform a numerical analysis showing the 
importance of our new terms. Finally, Section VI contains our conclusions. In the appendix 
supplementary material is provided: Appendix A contains the building blocks which are use- 
ful for the practical calculations and in Appendix B detailed analytical results are presented 
for the corrections to the matrix element {s^\0 2 \b}. For completeness the results of the 
order a s corrections are listed in Appendix C and intermediate results needed for the matrix 
element {s^\Oj\b) are given in Appendix D. In Appendix E the results are provided which 
are necessary to discuss the branching ratio BR(6 — > A s 7) S7 > Scut where E cut represents a 
cut-off on the photon energy. 

II. VIRTUAL CORRECTIONS TO b -> s 7 ASSOCIATED WITH 1 AND 2 

In this section we derive the (renormalized) order a 2 corrections to the matrix elements 
(s7|Oi|6) and {s^\0 2 \b). Thereby only the contributions proportional to the number of 
fermion flavors, rif, are taken into account. We show at the end of this section that the result 
for (s~f\Oi\b) can easily be obtained from the one for (s^\0 2 \b). Therefore, we concentrate 
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in the following on the calculation of the renormalized matrix element M2, 

M 2 = (s 7 |0 2 |6) , (5) 
which is conveniently written in the form 

M 2 = M 2 (0) + M 2 (1) + M 2 (2) . (6) 

The superscript counts the factors of a s . The leading term vanishes, i.e. = and 

the 0(a s ) calculation has been performed in [20]. In the following, we discuss the 0{a 2 a nf) 
term, In Subsection A we present the calculation and results of the dimensionally 

regularized three-loop diagrams, while Subsection B is devoted to the calculation of the 
counterterms. In Subsection C we combine the results of the three-loop results with the 
counterterms and derive the renormalized expression M 2 . 



A. Regularized three-loop corrections to (S7IO2I&) 

(2) 

The three-loop diagrams contributing to M 2 can be divided into four non-vanishing 

classes as shown in Figs. 1 and 2 1 . The sum of the diagrams in each class is gauge invari- 

(2) 

ant. The contributions to the matrix element M 2 of the individual classes are denoted by 
<L(1), <bL(2), <Le(3) and M 2 (2 b } arc (4) , where, e.g., M 2 (2 b } arc (l) is 

<L(1) = <L(1«) + <L(U) + <L(lc) • (7) 

For the practical calculation we essentially follow the techniques developed in [20]. To 
make the paper self-contained, we nevertheless present as an example the calculation of the 
diagram 2c in some detail. 

The amplitude ^2 bare (2c) is constructed with the help of the building blocks Ip and 
K^p,, shown in Fig. 9 in Appendix A. The analytic expression for Ip is given in Eq. (Al), 
while Kpp, is given in Eq. (A2) for an arbitrary mass m/ of the quark in the loop. This 
mass is retained in Kpp,, because it will be used as a regulator of infrared singularities in 
the calculation of (S7IO7I&). As (s 7 |02|6) is free of infrared singularities, we can put in this 
section mj = 0. Thus the parameter integral in Eq. (A2) can be expressed in terms of Euler 



In principle there are also diagrams in which the photon is emitted from the quark-loop insertion in the 
gluon propagator. However, these contributions vanish due to Furry's theorem. 
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la 





lb 





lc 




2a 



2b 



2c 



FIG. 1: Diagrams la-c and 2a-c associated with the operator O2. The photon is represented by 
a wavy line and is emitted from a down-type quark in all the diagrams. The virtual gluons are 
represented by curly lines. The sum of the first three graphs is denoted with ^2 bare (-0' whereas 
the sum of the second three diagrams is called Af^barc^)- 



T functions. Furthermore, only the gpp term has to be kept as the other building block Ip 
is transversal. The diagram 2c can be written as 

2j /as 

(47T) e V 7T 

( jd r 



M (2) 



2, bare 



(2c) 



eQ d C F l n f — 



(27T) d 

jf + i + m b 



-J- 



r(4 - 2e) 
■ft + / + m b 



(p' + r) 2 —ml + id (p + r) 2 — m 2 + i5 



-7 / ^M(p) 



(r 2 + i5) 1+ * 



dxx 



l-e, 



1 — X 



,1-e 



r — 



+ i5 



(8) 



, v x(l - x) 

where u{p) and u{p') are the Dirac spinors of the b and s quark, respectively, while e 
denotes the polarization vector of the photon. Cp and T are the eigenvalue of the quadratic 
Casimir operator and the index of the fundamental representation of the color gauge group, 
respectively, with the numerical values Cf = 4/3 and T = 1/2. The Euler constant 
7£ appears in Eq. (8), because we write the square of the renormalization scale in the 
form /i 2 exp(7 S )/(47r). The parameter 5 (with 5 > 0) in the denominators of the various 
propagators symbolizes the "e-prescription" . 




FIG. 2: Diagrams 3a-b and 4a-b associated with the operator 2 . The photon is represented by 
a wavy line and is emitted from an up-type quark in all the diagrams. The virtual gluons are 

(2) 

represented by curly lines. The sum of the first two graphs is denoted with M 2 bare(^)' whereas the 
sum of the second two diagrams is called M\ ^ are (4) . 

In a next step we denote the four different denominators with 

Di = (j/ + r) 2 - ml + i 5, 
D 2 = (p + r) 2 -m 2 b +i5, 

x{l — x) 
D 4 = r 2 + i5, 

and introduce a Feynman parametrization as follows: 

1 r(3 + 2e) r dudvdyw e y e ~ l 

D 1 D 2 DlDl +e ~ r(e)r(l + e) J {D lU + D 2 v + D 3 y + D A wf+ 2 ^ i " i 

with w — 1 — u — v — y. The integration variables (u, v and y) run in the simplex S defined 
through u, v, y > and u + v + y < 1. After the integration over r one simplifies the 
remaining integrals with the help of the substitutions 

u^{l-u')(l- l -^f] , v ^ 1^.(1 y^u'yi. (10) 
\ u J u 

The integration variable v' varies in the interval [1 — u', 1] whereas the other three variables 
x, y' and u' all vary in the interval [0, 1]. We tighten the notation by omitting the primes 
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and arrive at 

OA 1 ( as Yrnr r n r(^)r 2 (2-6)(i-6) 

JO JO Jo Jl-u 

(1 - yfu^uip') \J±- + A + 5^ ] «(p) , (11) 

where the Dirac matrices Pi, P2 an d -P3 are polynoms in the Feynman parameters and the 
expression C is given by 

C = m 2 b {l - u)v + -^—m\ - i5 . (12) 

We should mention at this point that the expression in Eq. (11) is infrared finite and is 
therefore regularized for e > 0. 

We use the same approach as in [20, 26, 27] and introduce Mellin-Barnes representations 
for the denominators C 1+3t and C 3t . In general the Mellin-Barnes representation of an 
expression of the form (K 2 — M 2 )~ x (with A > 0) reads 

(x?^=(Wrki/ ds (-^)' r( - s)r(A + s) - (13) 

where the integration path 7 runs parallel to the imaginary axis. It intersects the real axis 
somewhere between —A and 0. The Mellin-Barnes representation for C x , (A G {3e, 1 + 3e}) 
is implemented by identifying K 2 and M 2 as 

K 2 <-> m 2 (l — u)v , 

M 2 <-> r ^y m 2 c +t5. (14) 

x{l-x) 

The integration path 7 has to be chosen such that the parameter integrals exist for all values 
of s G 7. This means in our case that 7 has to intersect the real s-axis between — 3e and 
0. After interchanging the order of integration, the four Feynman parameter integrals can 
easily be expressed in terms of products of Euler T-functions. What remains to be done is 
the integration over 7 in the complex s-plane. We close the integration path in the right 
half-plane and use the residue theorem to perform this integral. The residues are located at 
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the following positions: 

s = 0,1,2,..., 
s = 1 — e, 2 — e, 3 — e, . . . , 
s = 1 - 2e, 2 - 2e, 3 - 2e, . . . , 
s = 1 — 3e, 2 — 3e, 3 — 3e, . . . , 

2 36, 2 36, 2 36, .... (l<-0 

The sum over the residues naturally leads to an expansion in the small parameter z = m 2 c /m% 
through the factor (rn 2 c /rnl) s in Eq. (13) (see also Eq. (14)). This expansion, however, is not 
a Taylor series because it also involves logarithms of z, which are generated by the expansion 
in e. The final result for Mfl &r . c (2c) can thus be written as 

<Lc(2c) = ^/^ fc ln i W, (16) 

k,l 

where the coefficients fk,i are independent of z. The power k is an (non-negative) integer 
multiple of \ and I G {0,1,2,3,4}. For a detailed explanation of the range of / we refer 
to [20]. 

In a similar way all other diagrams can be treated. The final result for the sum of the 
three-loop diagrams is given by 

<L = <L(1) + <L(2) + M 2 ( te(3) + M 2 ( te(4) , (17) 

where the analytical results for the individual terms of the r.h.s. are listed in Appendix B. 
We decided to include corrections up to 0(z 3 ) as the higher order terms lead to a negligible 
contribution for the physical value z ~ 0.1. 

B. Counterterm contributions to (S7IO2I&) 

In this section we work out the various counterterms of order a 2 rif which are needed 

(2) 

to derive the renormalized result M\ ■ There are counterterm contributions due to the 
renormalization of the strong coupling constant and due to the mixing of O2 into other 
operators. 

We first discuss the counterterms related to the renormalization of a s . As the leading 
term M2 ^ are is zero, only the renormalization of g s in the two-loop result generates 
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FIG. 3: Counterterm diagrams to O2 involving the operator O4. The crosses denote the possible 
places for photon emission. Note that the diagrams where the photon is emitted from the fermion- 
loop are zero due to Furry's theorem. 



a counterterm which can be written as 



5Z, 



(i)> n / 

9s 



25Z, 



(iW M (D 



9s 



2, bare ' 



a s TlfT 

7r 6e 



(18) 



M. 



(i) 



2, bare 



is the sum of the two-loop diagrams which has to be known including terms of 0(e). 
For this reason we extended the calculation of Ref. [20] to order e 1 . 

We now turn to the counterterms induced through the mixing of 2 with other operators. 
First, we consider the counterterms connected with the mixing of O2 into four-fermion 
operators. At order a s there are non- vanishing mixings into 0±, 4 and into the evanescent 
operator P u , denned in Appendix A of Ref. [12]. As only 4 has a non- vanishing matrix 
element for b — > 57 proportional to a s nf, the only counterterm of this type is 



M. 



(2) 



24,a 
,(1) 



SZ^M^ 
a s 1 



0Z, 24 — 7~ ) 

it be 
1 



Mi 1} = — I - — + 78 + 288 In ( — ) + 36m + 1159e - 1507r 2 e 



72 



81 



/i 



m b 



-312 In ( e - 576 ln z e + 258i7re - 144i7r In ^ ) e + C(e 2 ) 



T i C F Tn / Q d (s7|0 7 |6)t r ee , 

47T 



(19) 



where can ^ e f° un d m [12]. The Feynman diagrams contributing to i.e. to the 

corrections of 0(a s nf) to (s7|0 4 |6)t re e, are depicted in Fig. 3. They were computed following 
the strategy outlined in Section II A. 

At order a^nf, there are mixings of 2 into Oi, 4 and P n and again only 4 has a 
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matrix element of 0(af). Thus the only counterterm of this type reads 



SZ 



M (2) 

(2),n/ 
24 



5Z. 



(2),"/ w (0) 
24 iW 4 

18? ■ 



M, 



(0) 



^ — fi) 



e + 



12 



+ 2 In 2 



/i 



6 2 + 



C(e 3 ))c F Q d ( S7 |C 7 |6) tree . (20) 



In a second step we consider the counterterms connected with the mixing of O2 into the 
dipole operators 7 and 0$. One can easily see that only one counterterm of this type 
generates a contribution of 0(a^nf): O2 mixes at three-loop order into O7; in turn, from 
O7 the tree-level matrix element for b — > 57 is taken. The resulting counterterm therefore 
reads [12, 28] 



M 2 (2) = 5Z^' nf ( S1 \0 7 \b) tree , 
6Z%»< = {^) 2 C F Tn f 



Qu 

24 



81 



Qu 2Q C 



(21) 



where Q u = 
respectively. 



144 243 

2/3 and Q d = —1/3 are the charge factors of up- and down- type quarks, 



C. Renormalized result for (s7|C>2|b) 

Combining the three-loop result Mfl &re , calculated in Subsection A, with the various 
counterterm contributions discussed in Subsection B (see Eqs. (18), (19), (20), and (21)), 
we get an ultraviolet finite result. As mentioned earlier, the result is also free of infrared 
singularities. Inserting the numerical values for the color factors (Cf = 4/3, T = 1/2) and 
the electric charge factors (Q u = 2/3, Qd = —1/3), we get the following renormalized result 
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with 

,(2) 



- (23) 

Re (/f) = ( - 145 + ( 288 ~ 307r2 ~ 216 C(3) + 216L - 54vr 2 L + 18L 2 
+6L 3 ) z + 24ttV /2 + 6 (18 + 2vr 2 + 12L - 6tt 2 L + L 3 ) z 2 

- (9 + 14tt 2 - 182L + 126L 2 ) z 3 ) + 0(z 4 ) , (24) 

Im (^4 2) ) = ^ ( - 22 + (180 - 12vr 2 + 36L + 36L 2 ) z 

- (12tt 2 - 36L 2 ) z 2 + (112 - 48L) z 3 ) + £>(z 4 ) , (25) 



Re frf ) = ^ - 1^- - JL (H280 - 1520tt 2 - 171vr 4 - 5760C(3) 
V 2 / 6561 729 1215 v w 

+6840L - 1440tt 2 L - 2520C(3)L + 120L 2 + 100L 3 - 30L 4 ) z 

(43 _ 12 ln(2) - 3L) z 3/2 - -J- (11475 - 380tt 2 + 96vr 4 

+7200C(3) - 1110L - 1560tt 2 L + 1440C(3)L + 990L 2 + 260L 3 

_ 60 L 4 ) z 2 + ^1^5/2 _ _JL ( 6 2471 - 2424vr 2 - 33264C(3) 
; 243 2187 V SV ; 

-19494L - 504tt 2 L - 5184L 2 + 2160L 3 ) z 3 + 0(z 7/2 ) , (26) 
Im (rf } ) = ^ (495 - 12 (375 - 19tt 2 + 36C(3) + 84L + 48L 2 - 6L 3 ) z 
+6 (207 + 38tt 2 - 72C(3) - 126L - 78L 2 + 12L 3 ) z 2 
+8 (67 - 12vr 2 - 48L) z 3 ) + 0{z A ) , (27) 

where L = lnz. We note that in the derivation of this 0(a 2 nf) result, there was no need 
to renormalize the parameter in the corresponding 0(a\) expression. Therefore, the 
symbol (s7|07|6) tr ec can be interpreted to be (in and M^) 

g 

(s7|0 7 |6) trcc = m b —u(p')^u(p) , (28) 

where nib denotes the pole mass of the b quark. Concerning this point, the reader is also 
referred to Section III. 

We now turn to the renormalized matrix element M{ , associated with the operator 0\. 
Oi, defined in Eq. (4), can be written as 

1 = ^0 1 -^0 2 , (29) 
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(a) (b) (c) 

FIG. 4: Virtual (a), gluon-bremsstrahlung (b) and quark-pair radiation (c) graphs to the operator 
O7. In (b) and (c), the diagrams where the gluon is emitted from the s-quark are not shown. 

with 

0i = {Hi A) (4i^t) , (30) 

where a and (5 are color indices. It is easy to see that Oi has a vanishing matrix element 
for b — > S7. Therefore, one obtains 

M[ 2) = -^M 2 (2) . (31) 

III. REAL AND VIRTUAL CORRECTIONS TO (s^\O r \b) 

In this section we describe in some detail the steps needed for the calculation of the 
0(a 2 s rif) corrections to the matrix element (s^f\Or\b). Due to the presence of infrared singu- 
larities, the practical calculation proceeds in a slightly different way than for 0<i- As these 
singularities only get canceled at the level of the decay width when combining the virtual 
corrections shown in Fig. 4(a) with the gluon bremsstrahlung (Fig. 4(b)) and the quark-pair 
emission process (Fig. 4(c)), we first derive expressions for the O(alnf) corrections to these 
three contributions to the decay width. The corresponding expressions necessary to evaluate 
BR(_B — > Xg^E^Ecut are discussed in Appendix E. 

To fix the notation, we write the contribution from O7 to the decay width F(b — > X s ^) 

as 

F° 

1 w 



77 - 1 77 



1 + f w + f &w 



77 ~ 32tt 4 ^ 7 ' ' ^ ' 



The 0(a s ) correction, f 7 ^, can be extracted from Ref. [20], reading 

a s ( 32 16tt 2 64 (m b \ , 

r77 = 47 v -Y-^ + y ln U ,) • (33) 
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We further split ff 7 ),n/ in Eq. (32) as 

f (2),n f = f (2),(a) + f (2),(6) + f (2),(c) ^ ^ 

with obvious notation (Fig. 4). 

For the calculation of the three parts contributing to f ^' Uf we could in principle put 
rrif = m s = at the beginning of the calculation and use dimensional regularization for 
both infrared and ultraviolet singularities. We found it easier, however, to use the strange 
quark mass, m s , and the mass of the quark in the fermion bubble, rrif, as infrared regulators. 
For formulating the results, it is convenient to introduce the dimensionless quantities 

m 2 rn 2 f 

r = — 2 , /=-£. 35 
mi mi 

We now turn to the calculations of T^^ c \ f and f (in this order). 

Inspecting the explicit expressions for the quark-pair radiation process (cf. Fig. 4(c)), one 
finds that it can be worked out in our "massive" regularization scheme in d = 4 dimensions. 
Furthermore, one realizes that one can also put m s = 0, provided m/ is kept at a (small) 
fixed value. As a consequence, the quark-pair radiation process is completely regularized 
by the mass m/. The evaluation of this process is quite standard: in a first step the 
subprocess b — > s'-yg* is considered where g* represents a virtual gluon. Subsequently the 
other subprocess, describing the decay of g* into two fermions, is added. It is straightforward 
to perform the occurring phase space integrations where only the one over the gluon virtuality 
is non-trivial. However, in the limit mj — > also this one can be performed analytically. 
One arrives at the following result for the quark-pair emission process: 



fgMc) = /«s\ T_i2662 + 24vr 2 + 2592((3) + (144tt 2 - 5916) ln(/) 
V47T/ 243 

-9001n 2 (/)-721n 3 (/)] . (36) 



Due to the Kinoshita-Lee-Nauenberg theorem, it follows that the sum of the virtual and the 
gluon bremsstrahlung corrections also must be finite for d — > 4 and m s — > for fixed mj. 

We now turn to the gluon bremsstrahlung process. The diagram in Fig. 4(b) (combined 
with the one where the gluon is emitted from the s-quark) can be written as 



M<i>f = ^M«*>, (37) 
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where M^'^ denotes the lowest order matrix element for b — > s^g and 5Z^' nf reads 



sz. 



(!)>«/ 



.^?-(lH-241n ^+^ + 241n 2 
7r36\e y /i J \ jj, 



e + 0(t 



(38) 



Note that the 1/e pole is of ultraviolet origin; the infrared singularity is regulated by mj in 
this expression. In addition, there is a counterterm contribution due to the MS renormal- 
ization of the strong coupling constant of the form 



M. 



7,ct 



(39) 



with 



, 7 (i),n f _ a s n f T 

0/j 9s — 7. — 

7r 6e 



(40) 



Combining M.f^l with Mj^ b \ one obtains the renormalized matrix element M^' 



(b) 



M. 



J 9s 



+ 



M. 



(41) 



from which the 0(a 2 nf) contribution to the decay width is obtained in a straightforward 
way. One gets 



,(2),(6) 
77 



2 i 5zf}^ 



SZ. 



r 



(!),(&) 
77 



a s \ 2 CfTuj 
4ttJ 18 



^2 ln(/) + 4 In ^ ) + ln(/) ln(r) 



+21n( ln(r 



8tt 2 + 416 ln(/) - 32/r 2 ln(/) 



is ln 2 (/) + 832 In ( ^ - 647r 2 ln - 960 In 2 



-576 ln(/) In 



-36 In 



fi 

- 120 In 



41n(r)( rr 2 -181n(/) + 61n 2 (/) 
2 ( m b 



241n 2 (r) ln(/) + 2 In 



+ 721n(/)ln( ^ 
fi J \ fi 



(42) 



where f^'^ is the corresponding (normalized) decay width for b — > s^g in the 0(a s ) 
approximation. As in our regularization scheme the sum 5Z^' nf + 5zi 1] ' nf /2 is finite (in e), 
f$'^ is only needed up to terms of order e°, which simplified the calculation. 
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We now turn to the evaluation of the virtual corrections shown in Fig. 4(a) and also 
discuss the various counterterm contributions. For the diagram shown in this figure, we 
obtain 

r(2),(o) _ 1 



M, 



7, bare 



81 



%r (2 ln(r) - 1) + — ( 2 + 12 ln(r) - 6 ln(r) ln(/) 



-24 ln(r) In f^j - 3 ln 2 (r) + 6 ln(/) + 12 In (^j 
+ 1718 + 123rr 2 + 840 ln(/) + 36rr 2 ln(/) + 90 ln 2 (/) 



+ 181n 3 (/) - 144 In 



432 In 2 



-432 ln(/) In J + 18 ln ( r ) (^ 24 + * ~ 12 ln (/) + 3 ln (/) 
-48 In (y)+48 In 2 ^+241n(/)ln^ 



-:,Jbr(r) | 2 ln( f) - ( y ) : - 18hi :i (r) 



a 



|?r ) C F Tn f ( S1 \0 7 \b} tree . 



(43) 



We stress that this expression is derived in such a way that m s is understood to be sent 
to zero prior to rrif. This procedure is justified by the fact that for fixed rrif the sum of 
the virtual- and gluon bremsstrahlung contributions must be finite in the limit m s — > 0, as 
discussed above. 

The counterterm contribution M^ 2 ^ at 0(a^rif) has various sources. There is a contri- 
bution M^l^ due to the renormalization of g s in the 0(a s ) vertex diagram (i.e. like the 
one in Fig. 4(a), but without the fermion bubble), yielding 



iw 7,cti 



12 6 / / Tfih 

— ln(r) ln(r) I 4 — ln(r) — 4 ln I — 



+12 - ln(r) ( 48 + tt 2 - 48 ln ( '-^ ) - 2 J hi 



/' 

2 ( ™* 
/i 



+ ]2hr(/-) | ! - hi ( y 



ol \ ^ 

C F Tn f (sj\O r \b) tlee . 

471"/ 



21n 3 (r) 



Then, there is a counterterm contribution M^J^ 1 of the form 



M. 



(2),(a) 



SZ. 



(2),n/ 
2,6 



7,ct 2 



+ 



+ *z£ )>B ' + 5Zr 6 ' (2) ' n/ ( S 7|0 7 |6) tree . 



(44) 



(45) 
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Here, 5Z^ b ' nf and 5zf^ nf are the 0(a 2 nf) pieces of the on-shell wave function renormal- 
ization constants for the b and s quark, respectively, while the operator renormalization 



factor SZ 7 2 )' nf refers to the MS scheme. Note that the presence of the on-shell renormaliza- 
tion factor 5Zm^ h ' L) in Eq. (45) implies that in the lower order contributions the symbol 



7 on,(2),ny 



(s7|07|6)tree is understood to be the tree-level matrix element of 7 in which the running 
6-quark mass is replaced by the corresponding pole mass. The explicit form of the various 
SZ factors occurring in Eq. (45) can be seen in Appendix D. 

After combining Eqs. (43), (44) and (45) into the renormalized matrix element, the 
calculation of f is straightforward. We obtain 

-216 



p( 2 ),(a) _ fa s \ 2 C F Tn f 
177 " 



21n(/)+41n( — 
fi 



ln(/)ln(r) 



+2 In {^^j ln(r)^ + 7495 + 624tt 2 + 1086 ln(/) + 72tt 2 ln(/) 
+666 ln 2 (/) + 36 ln 3 (/) - 6336 In 



+ 6048 In' 



+2592 ln(/) In ( ) + 18 ln(r) ( vr 2 - 18 ln(/) + 6 ln 2 (/) 



-36 In | ^ 



m b 



// 



120 In 2 ^ + 721n(/)ln — 



108 br(r) | 1h(/H 2 hi ( ^ 



(46) 



We now combine virtual and gluon bremsstrahlung corrections given in Eqs. (46) and (42), 
respectively. We obtain (after putting T = 1/2 and C F — 4/3) 

V4tJ 243 

+721n 3 (/) - 576(9 + vr 2 ) In 



n(2),(a)+(6) 
77 



14990 + 1176tt 2 + (5916 - 144rr 2 ) ln(/) + 900 ln 2 (/) 



+ 3456 In 2 



m b 
fi 



(47) 



where the 1/e poles and the mass singularities associated with m s are canceled. 

When combining this result with the quark-pair emission process in Eq. (36), we obtain 
the final result 



,(2),n/ 



77 



477/ 



n f I 2tf ] In 



m b 
fi 



+ 2lf ] In 



m b 
fi 



+ 2r ( 7 



(2) 



(48) 
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FIG. 5: Graphs associated with virtual corrections to the operator Og- The crosses denote the 
possible places where the photon can be emitted. 



with 



,(2) _ 64 
h ~ 9 ' 
(2) 32 

4 



4 2) = -i(9+* 2 ). 



r 



S1 (97 + 50tt 2 + 108C(3)) . (49) 

The cancelation of the ln(/) terms is a strong check for the correctness of the individual 
pieces of the calculation. 

For later convenience we formally introduce an amplitude M 7 in such a way that its 
square reproduces the result of Eq. (48). Adopting the notation of Eq. (6) one gets 

M? = (g) 2 n f ( S1 \O r \b) tree (4 2) m 2 (^) + I? In (^) + rf) . (50) 

IV. VIRTUAL CORRECTIONS TO (s 7 |0 8 |6) 

We first discuss the two-loop diagrams depicted in Fig. 5, which contain the building 
block Kpp, (see Eq. (A2)). As these diagrams are free of infrared singularities, we put the 
masses of the quarks in the fermion loop as well as the strange quark mass m s to zero 
from the beginning. The calculation can be performed along the same lines as described 
in Section II A. However, due to the absence of m c , the actual evaluation of the diagrams 
turns out to be much simpler. The result can be cast into the form 

-4e 



KL = {^) 2 C F Tn f Q d ( S1 \0 7 \b} tTCC ^^ + 1(120-6^ + 18*70) (^) 



+530 - 28tt 2 - 180C(3) + 93iir 



(51) 
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The counterterm contribution of O(a^nf), denoted by Mg 2 ^, stems from the renormalization 
of g s and from the mixing of 8 into the operator 7 . We obtain 



Afg = SZ£» nf (si\0 7 \b) tICC + 25 M?> , 



r(2),n 



(52) 



with 



(2),n/ 
87 



(?r^i(?- 7 )' 



(l),n/ _ Os_rifT_ 

7T 6e 



3s 



M a (1) = 



fts i 

■T i o < 3dCj!'(s7|07|6)tree 



^ + 33-2^-24 ln(^ 



^G/tt + c | r2- Kr 2 -36C(3)-661n(^J + 4/r 2 In ^ 



2 / m & 



+241n z ( -^J + 12iTr - 12m In ^ ) ) f C(e 2 ) 



(53) 



SZ^' nf is obtained from [12, 28]. The sum of Mg 2 ^ are and Mg 2 ^ leads to the renormalized 
result (using T = 1/2, C F = 4/3 and Q d = -1/3) 



Mf ] = (^)' n f ( S1 \O r \b) 



tree 



£hn> r± +^ln ^ 



// 



,(2) 



(54) 



with 



,(2) 



64 
27' 



4 2) = f x (47 - 2vr 2 + 6m) , 



.(2) 



243 



(-314 + 16tt 2 + 72C(3) - 57m) . 



(55) 



V. NUMERICAL IMPACT OF THE 0(a 2 s n f ) CORRECTIONS 



It is well-known that the inclusive decay rate for B — > X s 7 is given by the corresponding 
6-quark decay rate Y{b — > X s ~f), up to power corrections of the form (AQ CD /m fe ) 2 [29] and 
(Aq CD /m c ) 2 [30] which numerically are well below 10%. 

As our new results are only a part of the complete NNLL contributions, we do not present 
a new prediction of the branching ratio in this paper. Instead, we only illustrate how the 
0{a 2 s nf) corrections to the matrix elements of the operators Oi, 2 , 7 and 8 modify the 
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NLL branching ratio for a given set of input parameters. For this purpose, we neglect power 
corrections and also electroweak terms. 

In a NLL calculation the inclusive quark-level transition b — > X s ^ involves the subpro- 
cesses b — > S7 (including virtual corrections) and b — > s^yg, i.e., the gluon bremsstrahlung 
process. We write the amplitude for the first subprocess similar as in Ref. [14]: 

AC 

^nll (6 ^ S7) = Vt * VtbD ^ {si lo 7 \b) tiee , (56) 

where the reduced amplitude _D NLL reads 

D NLL = Cf(fi) + ^ V^(pt) . (57) 

4.71 



The symbol V^(//), defined as 



i=l 



(58) 



incorporates the NLL corrections, and to the matrix elements. In Eq. (57), the first 
term on the r.h.s. is understood to be the Wilson coefficient C 7 s (fi) at NLL order, while 
the Wilson coefficients appearing in V^(fx) are understood to be taken at LL order. As in 
Ref. [14], we convert the running mass factor m b (fi), which appears in the definition of the 
operator O7 in Eq. (4), into the pole mass m&. This conversion is absorbed into the function 
V^(n) and consequently the symbol (s7|07|6) tr ce in Eq. (56) is the tree-level matrix element 
of the operator 7 , where the running mass factor mj(/i) is understood to be replaced by 
the pole mass mj,. The NLL virtual correction functions and in (58), taken from 
Ref. [20], are repeated for completeness in Appendix C. Note, that the quantity not only 
contains virtual corrections to the matrix element of 7 , which would be infrared singular. 
rj 1 ^ is constructed in such a way, that the (O7, O7) interference term generates the sum of 
virtual and bremsstrahlung corrections when formally calculating the branching ratio from 
^4 NLL (6 — > 57). For the details of this construction, we refer to Ref. [20]. Numerically, the 
square of this amplitude encodes the bulk of the decay width. The additional bremsstrahlung 
corrections, which are infrared finite for -E g i UO n — ► 0, are relatively small. Therefore, when 
considering terms of order 0(a^nf), we omit purely finite bremsstrahlung contributions. 

When improving the amplitude for the subprocess b — > 37 by including the terms of 
0(a 2 s 7if), the result can be written as 

AC 

A{b - si) = --^-V ts *V tb D( S1 \0 7 \b) tiee , (59) 
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where the reduced amplitude D is 



Air 



( 



An 



) 



2 



D = Cj 



(60) 



V^(/i), defined as 




(61) 



incorporates the 0(a 2 nf) corrections to the matrix elements calculated in the previous 
sections of this paper. The explicit Cj S (fx) term in Eq. (60) in principle stands for the NLL 
Wilson coefficient, supplemented by the rif dependent NNLL contributions. As the latter are 
not known yet, we take this Wilson coefficient at NLL precision in the numerical evaluations. 
The Wilson coefficients entering V^(fi) are in principle the LL coefficients, supplemented 
by the rif dependent NLL contributions. In practice, we decide to replace these Wilson 
coefficients by the respective complete NLL version. Finally, the Wilson coefficients entering 
V"( 2 )(/i) are the LL versions. Note, that the gluon bremsstrahlung and the quark- ant iquark 
emission processes associated with Oj are effectively transferred into and tf\ as 

described in Section III. As already mentioned above, the square of the so-defined amplitude 
incorporates the major part of the branching ratio. We therefore consider the additional 
finite bremsstrahlung corrections to the decay width only at the NLL level, i.e. we do not 
calculate the 0(a 2 rif) corrections to these contributions. 

As the square of the amplitude for b — > 57 (in the sense defined above) encodes the 
dominant part of the decay width, it is reasonable to compare the NLL result _D NLL in 
Eq. (57) with the corresponding (9(agn/)-improved result D in Eq. (60). In Fig. 6, the 
function D is plotted as a function of the renormalization scale \i. We note, as already 
discussed in the Introduction, that we use in the numerical evaluations the hypothesis of 
naive non-abelianization, which amounts to replacing rif by — 3/3 /2. Nevertheless, in the 
following we still write 0{a 2 s nf). The dash-dotted line shows the NLL approximation as 
defined in Eq. (57), while the solid curve shows the result after including the 0(a 2 rif) terms 
as discussed above. The dashed line shows the result with 0(a 2 rif) improvements, in which, 
however, the Wilson coefficients in (//) are taken in LL approximation. The three curves 
illustrate that the changes between the 0(a 2 rif) improved version (solid line) and the NLL 
prediction (dash-dotted line) are mainly due to the new 0{a 2 nf) corrections of the matrix 
elements calculated in the previous sections. 
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FIG. 6: The reduced amplitude D as a function of the renormalization scale [i where the plot 
on the right is an enlargement of the one on the left. The dash-dotted curve represents the NLL 
approximation and the solid curve includes the corrections of 0(a 2 nf). For comparison we also 
show the result where the Wilson coefficients in (cf. Eq. (58)) are inserted to LL precision 
only (dashed curve). 



From A(b — > 57) in Eq. (59) the decay width T(b — > 57) is easily obtained to be 



r ( fo - s 7) = t^j \V ts *V tb \ 2 a cm ml \D\ 



32tt 4 



(62) 



When giving numerical results for the NLL predictions, we only retain terms in \D\ 2 up 
to order a s , while for the improved version we retain terms up to 0{a 2 nf) in \D\ 2 and 
systematically dismiss higher order contributions. For completeness we should mention 
that a s (fj,) is evaluated using two-loop accuracy in the (5 function. We checked that the 
contribution of the three-loop term /3 2 is numerically small. 

To obtain the inclusive decay rate for b — > X s ^y, we have to take into account those 
terms which have not yet been absorbed into the virtual corrections. At NLL precision, 
these contributions consist of those gluon bremsstrahlung corrections which are finite when 
the gluon energy goes to zero; they have been calculated in Refs. [31, 32]. As the (0 8 ,0 8 ) 
contribution to T(b — > s^g) becomes infrared singular for soft photon energies, we introduce 
a photon energy cutoff E cut as in Ref. [12] and define the kinematical decay width 

r(6 - X sl ) Ej > Ecut . (63) 

At NLL the gluon bremsstrahlung contribution to this quantity can be written as 

r(6 - s ig ) Ei > Ezut = \V ts *V tb \ 2 a em m\ A , (64) 



327T 4 
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where A is of the form [12] 

TT 



A = (e-°.MMW+*MS))/W _ y + <tl £ Re [ C f (//) Cf (//) . (65) 

ij'=l;i<j 

The quantity 5 is defined through 

J Ecut = ! y(l-5)= J E max (l-(5). (66) 

In Eq. (65) we put C? ff = for % = 3, . . . , 6, as in the virtual contributions. We list the 
explicit expressions for the quantities fij(S) in Appendix C. 

We should repeat that the 0(a 2 nf) corrections are incorporated in the quantity D, de- 
fined in Eqs. (59) and (60). We stress that the absorbed gluon bremsstrahlung- and the 
quark-pair emission terms were obtained by integrating over the full range of the photon 
energy. Thus, since we decided to implement a photon energy cut as just described, the 
final expression for the kinematical decay width can be written as 

r(6 - X sl ) E ^ Ecut = ^\V ts *V tb \ 2 a cm m! (\D\ 2 + A) - r£>' B ' (b -> X sl ) Ej < Ecut , (67) 

where the expression for r^' n/ {b — > X s ^) Ei < Ecvit is derived in Appendix E. 
In a last step, the kinematical branching ratio is obtained as 

BR(6 - X sl ) Ej>Ecut = T{b - X '^^ BR SL , (68) 

1 SL 

where BRsl is the measured semileptonic branching ratio and the semileptonic decay width 
Tsl (supplemented by the 0(a 2 rif) terms [23]) is given by (z = m 2 /ml) 



-p G 2 F \V cb \ 2 ml 
FsL= 192.3 9(z) 



2a s (ji)h(z) f ajjj) \ ( fm c \ 1 h(z) ^ (m h 
37r g(z) \ 7r J °\ Xl3 \m b J 3g(z) R \ A* 



(69) 



where the phase space function g(z) and the 0(a s ) radiation function h(z) [33] read 
g(z) = 1 - 8z + 8z 3 - z 4 - 12z 2 ln(z) , 

Hz) = -(l-z 2 ) (|-^. + |^+,ln(,)(20 + 90.-^ 2 + ^ 3 

+ z 2 \n 2 (z) (36 + z 2 ) + (1 - z 2 ) - *±z + f , 2 ) In(l - z) 
- 4 (1 + 30 z 2 + z 4 ) \n(z) ln(l - z) - (1 + 16 z 2 + z 4 ) (6 Li 2 (z) - tt 2 ) 
tt 2 - 4 Li 2 (v^) + 4 Li 2 (-v/i) - 2 ln(z) In ( j— ^ 
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32^ 3/2 (l + 2) 



(70) 
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FIG. 7: The branching ratio as a function of the renormalization scale [i where the plot on the right 
is an enlargement of the one on the left. The dash-dotted curve represents the NLL approximation 
and the solid curve includes the corrections of 0{a 2 s rif). For illustration in the left plot the latter 
are also shown for the case where {M^) is set to zero which corresponds to short-dashed 

(long-dashed) curve. A photon energy cut of E cut = is used, which corresponds to 5 = 0.9. 

The function xpi m c/'mb), which encodes the 0(a 2 rif) terms 2 is given in the form of a 
plot in Ref. [23]. For m c /mb = 0.29, which is the default value in our paper, one finds 
X p(0.29) « 1.68. 

In Fig. 7 the kinematical branching ratio is shown for the choice E cnt = to&/20, or, 
equivalently, 5 = 0.9 [8] as a function of the renormalization scale //. The input parameters 
were chosen to be: m b = 4.8 GeV, m c /m b = 0.29, m t = 173.8 GeV, m w = 80.41 GeV, m z = 
91.187 GeV, a s (m z ) = 0.119, a cm = 1/137.036, \V ts *V th /V ch \ 2 = 0.95 and BR SL = 10.49%. 
The dash-dotted line shows the branching ratio BR(6 — > X s ^) in NLL precision. In this case 
the terms of O(a^/3o) are consistently omitted in the expression for Tsl in Eq. (69). The 
solid line shows the branching ratio where the C(o; 2 n/) (or the O(a 2 p )) improvements are 
included. 

One observes that for /j, xs 5.5 GeV the 0(a 2 s nj) corrections vanish and that they are 
negative (positive) for smaller (larger) values of //. In this context it is instructive to look 
at the decomposition of the result. For this reason we show in the left plot of Fig. 7 the 
0(a 2 nf) corrections where either and or is artificially set to zero which 

corresponds to the short-dashed and long-dashed curve, respectively. This illustrates that 

2 Note, that nf is replaced by — 3/3q/2. 
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FIG. 8: Dependence of the branching ratio on the photon energy cut, E cut = ^(1 — 5). The dash- 
dotted curve shows the NLL result, while the solid curve includes the 0(a 2 s nf) improvements. The 
renormalization scale is fi = 4.8 GeV. 

there is a large cancelation between the negative contribution from O-j and the one from 
0\ and 2 which is, of course, also present in the amplitude D. The effect of the a 2 nf 
corrections from the operator 8 is significantly smaller and at most of the order of 2% in 
the considered interval for fj,. 

Fig. 7 furthermore illustrates that the /i dependence of the 0(a 2 rif) improved prediction 
for the branching ratio is somewhat natter than in the NLL case if we restrict ourselves to 
fi > 4 GeV. This is a welcome feature of our result, however, in general we cannot expect 
to reduce the /i dependence as the solid curve only represents a part of the 0(a 2 ) result. 
Indeed, we obtain a stronger /x-dependence in the region below 4 GeV. 

In Fig. 8 we show the dependence of the kinematical branching ratio on the photon 
energy cut. The dash-dotted line shows the NLL result, while the solid curve includes the 
order a 2 nf improvements. We should mention at this point that we did not include any 
non-perturbative effects in the photon energy spectrum. The main purpose of this figure is 
to illustrate how the order a 2 nj contributions modify the NLL result. 
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VI. CONCLUSIONS 



In this paper a first step towards a complete NNLL calculation is undertaken and radiative 
corrections to the matrix elements of the operators 1 , 2 , 7 and 8 are computed. More 
precisely, we consider the contributions of order a 2 rif which are induced by a massless 
quark loop. It is expected that these corrections, after replacing rif by — 3/?o/2, may give 
an important contribution to the full order a 2 corrections. Furthermore, motivated by the 
NLL analysis, we expect that the 0(a 2 nf) corrections to the matrix elements numerically 
dominate the ones of the same order to the Wilson coefficient functions and to the anomalous 
dimension matrix. 

In practice our calculation requires the evaluation of two- and three-loop diagrams in the 
case of O7, 0$ and Oi, O2, respectively. Furthermore, in order to obtain an infrared finite 
result in the case of O7, also the contributions from the gluon bremsstrahlung and from the 
quark-pair emission process are taken into account which requires the evaluation of three- 
and four-particle phase space integrals, respectively. All calculations are performed analyt- 
ically where an expansion in m c /m b is applied to the three-loop diagrams. For practical 
purposes this expansion is equivalent to the exact result. 

As far as the numerical impact of our result is concerned, we observe a striking cancelation 
among the individual contributions at order a 2 nf. When using a photon energy cut of 
-^cut — 771&/20, the 0(a 2 rif) terms reduce (after replacing nf by — 3/3 /2) the branching 
ratio by —0.98% for /j, = m b = 4.8 GeV and lead to corrections of —3.9% and +3.4% for 
/j, = 3.0 GeV and pL = 9.6 GeV, respectively. 
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FIG. 9: The building blocks Ip and Kip, which are used in the calculation of the Feynman diagrams. 
The curly lines represent virtual gluons, whereas the letters b, c and s stand for the corresponding 
quark (/ stands for a generic quark of mass rrif). Note that the external gluons are not amputated 
in the case of K^,. 





Ja/3 

FIG. 10: The building block J a p used in the calculation of the Feynman diagrams involving 0\ 
and O2. The curly and wavy lines represent off-shell gluons and on-shell photons, respectively. 

APPENDIX A: BUILDING BLOCKS 



The three-loop diagrams involving 0\ and O2 as well as the two-loop graphs involving 
O-j and 8 can be calculated by using one or more of the building blocks Ip, J a p and Kip, 
to be discussed in this appendix. The corresponding diagrams are shown in Figs. 9 and 10 
where the color indices are suppressed. 

The calculation of Ip is straightforward and yields 

A 



^r( £ )/i*ew(l- £ ) e *™ {rpf-r^p) L — 



f 

Jo 



dx [x(l — x)] 



l-e 



x(l — x) 



iS 



(Al) 
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where r is the momentum of the virtual gluon emitted from the c-quark loop. In the three- 
loop diagrams shown in Fig. 1 (cf. Section II A), the free index (3 will be contracted with 
the corresponding index of the dressed gluon propagator Kim. 

It is also quite simple to obtain the building block K^, (i.e., the dressed gluon propagator) 
which can be cast into the form 

2 1 _ r P r p' rl 

Rf = -|^rr(e)e^VV £ - ^' 2 T~ \ &xx{l-x)(x{l-x)r 2 -m) + i5)-\ 

(A2) 

where rrif denotes the mass of the quarks and T —\. Note that this expression is indepen- 
dent of the gauge parameter £ which enters the free gluon propagators in the construction 
of Kpp, when working in an arbitrary gauge. 

The building block J a p is somewhat more involved. Adopting the notation of Ref . [34] , 
it reads (for an on-shell photon) 

^ 9s Qu 



Ja/3 — 



16 7T 2 



E(a, [3, r)Ai 5 + E(a, 0, q)Ai 6 - E{/3, r, q)— Ai 23 

q ■ r 



-E(a,r, q)^At 25 - E(a,r, q)^ Ai 26 



L \ , (A3) 



q ■ r q ■ r 

where q and r denote the momenta of the on-shell photon and the off-shell gluon, respectively. 
When inserted into the full diagrams in Fig. 2, the indices a and f3 will be contracted with the 
polarization vector e of the photon and with the dressed gluon propagator K^, , respectively. 
The matrix E(a,(3,r) is defined as 



E(a,(3,r) = I( 7a7/J /- TV*,), (A4) 
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and the dimensionally regularized quantities Aik occurring in Eq. (A3) read 

Ai 5 = AB + dxdy [4(q ■ r) xy (1 - x)e + r 2 x (1 - x)(l - 2x)e 
Js 

+(l-3x)C] C' 1 ' 6 , 

Ai e = 4 B + J dx dy [—4(5 • r) x y (1 — y)e — r 2 x{2 — 2x + 2xy — y)e 
-(l-3y)C] C- 1 " 8 , 



Ai 2 3 = — Ai 2 6 = 8 -B + (g • r) dxdyxyeC 



l.i u /y ,< /y i '.. ■ 1 € 

5 



Ai 25 = -8 J B + (g-r) dxdyx (1 - x) eC _1_£ , (A5) 

is 

where 5+ = (1 + e)r(e) e Tise /x 2e and C is given by 

C = m 2 — 2xy(q ■ r) — r 2 x (1 — x) — iS. 

The integration over the Feynman parameters x and y is restricted to the simplex S, i.e. 
y £ [0, 1 — x], x G [0, 1]. Due to Ward identities, the quantities Ai k are not independent of 
one another. Namely, 

(fJap = and r^Jaf] = 
imply that Ai 5 and Ai Q can be expressed as 

2 

Ai 5 = Ai 23 , Ai 6 = — Ai 25 + Ai 26 . (A6) 

5 • r 

APPENDIX B: REGULARIZED THREE-LOOP RESULTS FOR (s 7 |0 2 |6) 

In Section II A we explained in some detail the calculation of the virtual three-loop 
corrections to (s^|0 2 |&). Here we give the results for the four gauge-invariant sets of graphs 
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depicted in Figs. 1 and 2. The results read, using z = ml/ml and L = \n(z): 

<Le(l) = 



<L(2) = < 7 



1 29 1 1 

+ - (5 + 2L) z + - (5 - 2L + 2L 2 - 2n 2 ) 

81e 243 6 V ; 6 V 



2^ z 2 



\ 77T 

+ — (17 + 30L - 18L 2 + 18tt 2 ) z 3 - — (l - 9z + 9z 2 - 18Lz 2 



-10z 3 + 12Lz 3 ) 



-6c 



+ 



1063 19tt 2 
1458 + ~324 



+ -!- (61 + AL - 9L 2 - IOtt 2 ) z + (79 - 22L + 28L 2 - 8L 3 - 9tt 2 
18 18 

-IALtt 2 - 12C(3)) z 2 + (63 - 27L - 36L 2 + 24L 3 - 59tt 2 + 42Ltt 2 + 
36C(3)) z 3 (58 - AAlz - 9 (23 + 38L - 6L 2 - 12tt 2 ) z 2 

-12 (4 + 3L + 3L 2 + 6tt 2 ) z 3 ^ + 0(z A ) I C F Tn f Q d { S1 \0 7 \b) 



tree • 



(Bl) 



324 



(157 -6L- UAL 2 - 60tt 2 ) 



2^ z s 



-6c 



H5I + i^ + l (96 _ 17 ^_ 126C( 3)) 

1458 72 54 v sv ;; 



+ ^-(40-18L-721n(2))^ 3 / 2 

— I 

+ — (213 + 102L - 40L 2 + 8L 3 + 34vr 2 + 96C(3)) 



20tt 2 



,5/2 



324 



(2799 - 995L - 198L 2 + 192L 3 - IOtt 2 - 60Ltt 2 - 936C(3)) 



+0{z 7 ' 2 ) 1 (^) 2 C F Tn f Q d (s^\0 7 \b) tiee , 



(B2) 
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I J_ + l JL _ 1 (i 8 + 24L + 3L 2 + 2L 3 - 3vr 2 - 6Lvr 2 - 24C(3)) z 
e 36e 432 36 V SV 

■— (15 + 6L- 6L 2 + 2L 3 + 6tt 2 - 6Ln 2 - 24C(3)) ^ 2 

oD 

-— (17 - 12L) ^ 3 + — (3 - 24^ - 6Lz - 6L 2 z + 2tt 2 ^ - 6z 2 + 12Lz 2 
36 36 



-6LV + 2ttV - 12z 3 ) 



-6e 



6029 17tt 2 



2592 144 1080 



(7200 + 6240L - 120L 2 + 220L 3 - 105L 4 



-2040tt 2 - 1200Ltt 2 + 90LV + lllvr 4 - 4440C(3) + 1440LC(3)) z 

~9i^n ( 15135 ~ 5790L ~ 1050L2 + 980L3 ~ 210L4 ~ 307r2 ~ 780Ln 2 
+180LV + 222tt 4 - 4560C(3) + 2880LC(3)) z 2 + -|- (3 - 2L + 72tt 2 ) z 3 

i — 



ITT 



432 



(411 - 4 (786 + 192L + 93L 2 - 24L 3 - 49tt 2 - 12Ltt 2 - 72((3)) z 



+2 (309 + 102L - 186L 2 + 48L 3 - IOtt 2 + 24Lvr 2 + 144C(3)) 



+8 (75 - 54L) i 



0{z A 



7Y 



C F Tn f Q u (sj\0 7 \b) 

tree • 



(B3) 



1 
e 
1 



ik + I-s( 12+6L - L3 -' 2 - 3 ^ 2 - 12 « 3 ») 2 



(6 - 6L + 3L 2 - L 3 + 2tt 2 + 24((3)) 

6e 



36 

81L 2 - 27tt 2 ) z 3 
2839 13tt 2 



1 

324 



(27 + 108L 



1 



(9480 + 2040L + 180L 2 - 340L 3 + 105L 4 
2592 144 2160 v 

+260tt 2 - 720Ltt 2 + 30LV - 439tt 4 - 3360C(3) - 8640LC(3)) z 



8tt 



z 3/2 + — (29895 - 6270L - 1410L 2 + 740L 3 - 210L 4 + 920tt 2 
-480Ltt 2 + 120L V - 52tt 4 - 16320C(3) + 4320LC(3)) z 2 + ^~^' 2 



216 

+0{z 7 ' 2 



(1358 - 477L - 99L 2 + 90L 3 + 63tt 2 - 18Ltt 2 - 432C(3)) 



a \ 2 

C F Tn f Q u ( S1 \0 7 \b) tree . 

71 / 



(B4) 
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In these expressions, ( denotes the Riemann ( function with the value £(3) 1.2020569. 
Q u = 2/3 and Qd = —1/3 are the electric charge factors of the up- and down- type quarks, 
respectively, while C F = 4/3 and T = 1/2 are color factors. 



APPENDIX C: CORRECTION FUNCTIONS NEEDED FOR THE NLL RESULT 

The renormalization scale independent parts of the virtual corrections in NLL order 
precision, encoded in the functions rf\ appearing in Eq. (58), read 

r« = — {-833 + 1447TV/ 2 
2 243 1 

+ [1728 - 180tt 2 - 1296 C(3) + (1296 - 324tt 2 )L + 108L 2 + 36L 3 ] z 

+ [648 + 72tt 2 + (432 - 216tt 2 )L + 36L 3 ] z 2 

+ [-54 - 84tt 2 + 1092L - 756L 2 ] z 3 } 

+ 1^! |_5 + [45 - 3tt 2 + 9L + 9L 2 ] z + [-3vr 2 + 9L 2 ] z 2 + [28 - 12L] z 3 } 
81 

+ 0(z 7 / 2 ), 
(l) 32 8 2 

r7 = ' 

= -^(-33 + 27r 2 -6m), (CI) 

where z is defined as z = m 2 /m 2 and the symbol L denotes L = \n(z). The quantities l*f~\ 
appearing in Eq. (58), read 

/(i) _ 1 /(i) /(i) _ 416 /(i) _ 8 /(i) _ 32 (r9 \ 
l i — g f 2 > L 2 — ~^7~ ) 3 ' ~9~ ' 

Notice that r<p, r±~\ and r$~\ as well as l£\ lip and l^p are not needed in the 
approximation Cf ff (//) = (i — 3, 4, 5, 6). 

The functions needed for Eq. (65) are taken from Ref. [12] and are listed here for com- 
pleteness. Note that fn(8) differs from the one given in Ref. [12] in order to be compatible 
with our r 7 given in Eq. (CI) 3 . 

MS) = ±MS) , MS) = -\h 2 {5) , MS) = -IMS) , MS) = -iMS), 



3 The additional, (5-independent addend appearing in our fn{$) is such that 777(1) vanishes: the contribu- 
tion of /77(A) at S = 1 is already absorbed into our r-j. 
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/22(a) 

MS) 

fM 
/77(a) 

/tbO&) 

/88(5) 



I62 

^7 



(l-5)/ 2 

5 / dt(l-zt) 



.(!-*)/* 



G(t) 1 


to 

h 


t 2 


7(1 



.1/* 



(1 - ztf 



'(l-S)/z 



G(t) 1 
t 2 



5 / eft Re G(t) + - + / dt{\- zt)Re G(t) + - 

Jo V 2/ J(i-s)/z V 2 



-3/27^), 

^ + ^-^3 + l^_ 4)ln(5) _| 



vr 2 „,,,., 9. 1 „ 2 1 o 

2 1 5 3 



Li 2 (l-*)- y -51n(5) + i 5- i 5 ^ 



1 ' -21n(^ 
27 1 \m s 



5 2 + 2<5 + 41n(l-5) +4Li 2 (l-5)- 



2tt 2 



-5(2 + 5) \n(5) + 8 ln(l - 5) + 75 + 35 2 -^ 3 } 



where the function G(t) is defined through 



(C3) 



-2 arctan 2 { x hr-A for t < 4 

G(t) = <( VV 4 "' 7 (C4) 

+ 2 In 2 {\{V~t + Vt^I)) -2mln(§( v / t + v 7 ^ 7 ^)) for t > 4. 

The functions associated with the operators O3 — Oq are not needed in our approximation. 
Note that in the numerics we set m s equal to zero in all terms except f^(S), where a value 
of rub/iris = 50 is chosen. 

APPENDIX D: 0(a*n f ) CONTRIBUTIONS TO VARIOUS Z FACTORS 



In this appendix we give the results for the 0(a 2 nj) contributions for various Z factors 
entering the calculation of the counterterm M^f^ in Eq. (45) in Section III. For the 
meaning of the various terms, see the text after Eq. (45). The (9(a 2 n/) contributions to the 
relevant Z factors read 



<" = (7) ! ^(t('-«- ! 4^ 



[I 

+30tt 2 + 96 ln(/) + 72 ln 2 (/) - 264 In 

V A 1 

+2881n(/)ln(^) +432 In 2 ^ ' 



443 



(Dl) 
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sz. 



(2),n/ 



2,s 



5z: 



on,(2),n f 



SZ. 



(2),n/ 



77 



-44 ln(/) + 12 ln 2 (/) + 24 ln(/) ln(r) - 88 In 



- 5 + 2rr 2 
m b 



+961n(/)m^J + 48 ln(r) In j +1441n : ^ 

«- A 2 <vr n/ f 71 + g7r2 _ 104 ln (mA + 48 ln2 ^ 



/i 

2 / m 6 



(f) 



96 



+ - 



10 12 



2 y ' 



with r = m 2 /m 2 and / = m^/m 2 ,. 



(D2) 



(D3) 
(D4) 



APPENDIX E: IMPLEMENTING THE PHOTON ENERGY CUT-OFF IN THE 
0(a s n f ) TERMS 

In this appendix we provide the formulas which are needed to calculate the 0(a^nf) piece 
of the kinematical branching ratio BR(6 — > X s j) Ei > Ecut , where E CVLt represents a cut-off on 
the photon energy. As can be seen from the structure of Eq. (67), this amounts to calculate 
rfj' nf {b — > X s j) E ^< Ecut , which is contained in the quantity D of Eq. (67). 

Note that only the gluon bremsstrahlung- and the quark-pair emission processes enter the 
calculation for vfj' ni [b — > A s 7)£ 7 <_E cut as the photon energy in the virtual contributions is 
concentrated at m b /2. The 0(a 2 s nf) contribution to T^' nf (b — > X s ^) Ei < Ecxxt can be written 
in the form 



vT f {b^X sl ) E ^ t = Y^ 



f 77 \Ej < E cut ) + t[j'^ \Ej < E cut/ 



(El) 



where (b) and (c) denote the gluon bremsstrahlung- and the quark-pair emission process, 
respectively, and T® 7 is given in Eq. (32). Like in Section III we use a regulator mass rrif for 
the secondary quark- ant iquark pair which means that Eq. (El) can be calculated in d — 4 
dimensions and with m s = 0. 

The calculation for the gluon bremsstrahlung piece ff^ <yb \E 1 < E cut ) is straightforward. 
Adopting the notation 

(E2) 



E cut = ^(l-5)=E max (l-5), 
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the result reads 

t^ b \ El <E cut ) = (^) 2 ^^(31-305-35 2 + 25 3 + 211n(5) 

+125 ln(5) - 35 2 ln(5) + 6 ln 2 (5)) ^ln(/) + 2 In (^j ^ , (E3) 

with / = rrif/ml. 

The calculation for t^'^ c \E 1 < E cut ) is somewhat more involved but still can be per- 
formed analytically, yielding 

f(f c) (E 7 <E cut ) = (^) 2 ^^^- 147- 97r 2 + 48C(3)-48Li 3 (5) + 54Li 2 (5) 

- ln(5) (85 - 4tt 2 - 54 ln(l - 5) - 24Li 2 (5)) + 13 ln 2 (5) + 12 ln 3 (5) 
+5(l60 + 47r 2 -24Li 2 (5) -241n(l-5) ln(5) - 94 ln(5) + 36 ln 2 (5)) 

+5 2 ( 1 - 7T 2 + 6Li 2 (5) + 6 ln(l - 5) ln(5) + 19 ln(5) - 9 ln 2 (5) 



-5 6 \1A - 4 HS ) - 2 ln(/) ^31 - 305 - 35 2 + 25 3 + 21 ln(5) 

+ 125 ln(5) - 35 2 ln(5) + 6 ln 2 (5) ] ) . (E4) 



Note that the sum of f ^'^(E^ < E cut ) and Y^^ c \e i < E cut ) is finite in the limit m/ — > 0. 
This completes the calculation of rf^' n/ (b — > X s ~f)E^<E cut , defined in Eq. (El). 

We note that differentiating vf^ ,ni (b — > X s ^) Ej <E cut with respect to the photon energy 
cut _E cut generates the corresponding term of order a 2 n/ to the photon energy spectrum. 
The result we obtain is in complete agreement with Eq. (9) of Ref. [25], where 0(a 2 nf) 
corrections to the photon energy spectrum were calculated. 
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